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ABSTRACT 

Context. Optical emission is detected from filaments around the central galaxies of clusters of galaxies. These filaments 
have lengths of tens of kiloparsecs. The emission is possibly due to heating caused by the dissipation of mechanical 
energy and by cosmic ray induced ionisation. CO millimeter and submillimeter line emissions as well as H2 infrared 
emission originating in such filaments surrounding NGC 1275, the central galaxy of the Perseus cluster, have been 
detected. 

Aims. Our aim is to identify those molecular species, other than CO, that may emit detectable millimeter and sub- 
millimeter line features arising in these filaments, and to determine which of those species will produce emissions that 
might serve as diagnostics of the dissipation and cosmic ray induced ionisation. 

Methods. The time-dependent UCL photon-dominated region modelling code was used in the construction of steady- 
state models of molecular filamentary emission regions at appropriate pressures, for a range of dissipation and cosmic 
ray induced ionisation rates and incident radiation fields. 

Results. HCO^ and C2H emissions will potentially provide information about the cosmic ray induced ionisation rates 
in the filaments. HCN and, in particular, CN are species with millimeter and submillimeter lines that remain abundant 
in the warmest regions containing molecules. 

Conclusions. Detections of the galaxy cluster filaments in HCO"'', C2H, and CN emissions and further detections of them 
in HCN emissions would provide significant constraints on the dissipation and cosmic ray induced ionisation rates. 

Key words. Astrochemistry - ISM:abundances - galaxies: intergalactic medium - Galaxies: clusters: individuahPerseus 
- galaxies: individual: NGC 1275 - Galaxies: cooling flows 

1. Introduction countered in attempts to match the observed spectra with 

f — [ , r ^ the results of models in which photoionisation by stellar ra- 

I , The central galaxies of many clusters of galaxies are sur- diati on is assumed to be the dominant heating mechanism 

g . rounded by filaments observed m optical emission lines to (^ [Johnstone et al1[2007l) . 

, ' extend up to many tens 01 k iloparsecs (e. g. ICrawiord et al.l 

. . : fT999t IConseliceetaD[200l . The filaments in NGC 1275, Infrared H2 emission lines trace the optical emission 

> ■ the central galaxy of the P erseus cluster, are tvpica lly of "f the filaments (jJaffe et al. | |2005D. [ Salome et al. | ([20061, 

kp ' the order of 100 pc thick. 'McN amara ef^ (Il996h sug- l2QQ8aD obtained CO(l-O), and C0(2-l) single dish data on 

^ , gested that the filaments originate in interactions between NGC 1275 and its surrounding filaments and they observed 

C^ ' the hot intracluster medium and the relativistic plasma as- HCN(3-2) at the position of 3C84, the radio source located 

■ - - ■ sociated with the AGN-driv en radio feat ures found in and ^t the centre of the galaxy NGC 1275. These two sets of ob- 

around the central galaxies. iPooe et al.l (l2008b) modelled servations estabhshedthat both regions contain cold molec- 

the filaments of NGC 1275 as material ablated from clouds ulargas.l Lim et al. | (|2008D ; [ Salome et al. | (|2008bD;| Ho et al. | 

and dragged by a more diffuse outflo w moving away from (2009|) have obtained mterferometric observations of the 

NGC 1275 at more than 10^ km s-\ iPone et al.l toM^ C0(2-l) emission of some of the NGC 1275 filaments. 

suggested that dissipation occurring during the momen- The required heating rates per particle due to dis- 

tum transfer from the diffuse flow to the filaments might sipation and/or cos mic ray induced ionisation found by 

po wer the filamentary emiss ion, an idea related to that iFerland et al.l (J2009l ) are several orders of magnitude higher 

of ICrawford fc FabianI (119921 ) that turbule nt mixing lay- than those in the nearest Milky Way Giant Molecular 

ers pl a .y a r ole in generating the emission. iFerland et al.l Clo uds, and the require d thermal pressures identified 

(|2008l I2OO9D have constructed detailed models of the op- by 'Ferland et al." ('2009') (see also 'Sand ers et all 12004 

tical and infrared line emission of filamentary gas heated Sa nders fcJFabian 2007) are about an order of magnitude 

by similar dissipation and by cosmic ray induced ionisa- higher than those of the envelopes of dense cores associ- 

tion. They were motivated in part by the difficulties en- ated with low mass star formation. Though the heating 
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and ionisation rates in the cold molecular regions of the 
NGC 1275 filaments may be lower than those in the cor- 
responding optical emission regions, they still may greatly 
exceed those relevant to the nearby Milky Way dense cores. 
If these rates are relatively high, the chemical composition 
in the NGC 1275 filaments should differ from that in a typ- 
ical nearby dense core and may provide a diagnostic of the 
conditions in the filaments. 

Hence, we have examined the chemistry of molecular 
regions subject to heating rates per volume and cosmic ray 
ionisation rates per particle in ranges that extend up to 
those approaching those that, according to iFerland et al.l 
(|2009t ). obtain in the optical emission line regions of the 
NGC 1275 filaments. Our purpose is to identify those 
molecular species that may be detectable and to determine 
which of those species may act as diagnostics of the dissipa- 
tion heating and of the cosmic ray intensity in the filaments. 
We shall consider a range of parameters up to those indi- 
cated by the observations; however, it is not our purpose 
here to model the filaments in NGC 1275. 

Section [2] contains details of the physical and chemical 
models that we adopted to describe the molecular compo- 
nent of the filaments. In section [3] we present our results, 
and section 2] ends the paper with a discussion and some 
conclusions. 

2. The Models 

We compute the filament chemistry by assuming that the 
filament may be represented by a steady-state photon- 
dominated region (a PDR) with constant thermal pres- 
sure. To solve for the chemistry in the filament we em- 
pl oy the UC L_PDR time-depende nt code, as implemented 
bv lBeU et al . (2005, 2007) and Bav et et al.l (|2009) run for a 
sufficiently long time to reach steady state. This code oper- 
ates in one space dimension and computes self-consistently 
the chemistry and the temperature as functions of depth 
and time within a semi-infinite slab, taking account of a 
wide range of heating and cooling processes. In the present 
work, the code is used to determine the chemical and ther- 
mal properties at all depths up to a maximum of about 20 
visual magnitudes. For the present work, we have updated 
the UCL_PDR code to include additional radiative cooling 
due to rotational transitions of ^^CO, C^^O , CS, OH and 
H2O. The escape probability formalism of Ide Jong et al.l 
(|l980t ) is used to determine non-LTE level populations and 
resulting line intensities at each depth- and time-step, in 
the same manner as for the existing coolants in the code. 
Collisional rates with H2 are taken fro m the Leiden Atomi c 
and Molecular Database (LAMDA; ISchoier et alll2005[) . 
The inclusion of these coolants allows the thermal balance 
to be more accurately determined in warm dense gas at 
high extinction. The treatment of H2 self-shielding has also 
been upd ated to use the re sults of detailed calculations per- 
formed bv lLee et al.l(jl996|) . The chemical network links 131 
species in over 1700 gas-phase reactions; only H2 is formed 
by surface chemistry; freeze-out of species on to grain sur- 
faces is excluded. The UCL_PDR code has be en validated 
again st all other commonly used PDR codes (jRoUig et al.l 
|2007() . The main elemental abundances relative to hydro- 
gen that have been adopted allow for some depletion onto 
dust, and are C: 1.4 x lO""*; O: 3.2 x 10""; N: 6.5 x 10"^; S: 
1.4 X 10-6; Mg: 5.1 x 10"*^; and He: 7.5 x 10"^ The adopted 
gas: dust mass ratio is 100. 



The code requires the adoption of several important pa- 
rameters that define the local physical conditions. These are 
(i) the gas number density, or equivalently (since the tem- 
perature is determined by the code) the thermal pressure; 
(ii) the FUV radiation flux impinging on the slab bound- 
ary; and (iii) the cosmic ray induced ionisation rate, here 
assumed to be independent of space and time. Most impor- 
tantly, in this computation there is also an additional term: 
(iv) an assumed additional hea t ing so urce in the filaments, 
as indicated by iFerland et al.l (|2009D . We discuss each of 
these parameters in turn. 

(i) In most cases the thermal pressure was taken to be 
-P = ^b(4 X 10^ cm"'^ K), where ks is Boltzmann's con- 
stant. This value of P is c lose t o the the average values 
measured bv I Sanders et al.l ()2004D for fila ments showing op- 
tical line e mission in NGC 1275 (see also lSanders fc FabianI 
I2OO7D . The lSanders et al.] ()2004D measurements show a wide 
scatter, and the average shows an upward gradient in pres- 
sure with decreasing cluster radius. Therefore, in some com- 
putations we have set the pressure to be one order of mag- 
nitude larger than stated above. 

(ii) T he incide n t FUV radiation field was taken to be / 
times the iHabind (|1968f ) radiation field mean intensity in 
the interstellar medium of the Milky Way galaxy. Values 
of / of 10 and 100 were adopted, on the assumption that 
the central galaxy in the cluster is, like NGC 1275, a bright 
galaxy with a central AGN. Of course, obscuration of stel- 
lar and AGN radiation affects the intensity impinging on 
the filaments, and the nature and strength of the incident 
radiation field are uncertain. 

(iii) The cosmic ray induced ionisation rate, ^, was as- 
signed values in the range of 5 x 10~^^ s~^ to 5 x 10~^^ 
s~^; th e higher values reflect the suggestion of lFerland et al.l 
(|2009D for filaments in NGC 1275. Values around 5 x 10"!^ 
s"""^ are often assumed for clouds within several hundred 
parsecs of the Sun and even in some studies of extragalac- 
tic regions. 

(iv) A special feature and most important modification 
of the present calculations is the inclusion of an unspec- 
ified heat source, in addition to the familiar main heat 
sources of photo-absorption and cosmic ray induced ionisa- 
tion. The add i tional heating mechanism, is here - following 
IFerland et al.l (|2009 f) - assumed to be dissipation. We in- 
corporate it by adding a heat source term, H, to the equa- 
tion of thermal balance. It is time-independent and uni- 
form everywhere, and specifies the thermal energy added 
per unit time per u nit volume by processes no t otherwise 
included in the code. IFerland et al.l (|2008l |2009D found that 



H must be within an order of magnitude of 10~^^ erg cm^'^ 
s~^ if dissipation powers t he observed optical emission re- 
gions. l^rlan^^F^Ll (|2009l) chose to make their heat source 
term proportional to gas number density. However, since 
the physics of the dissipation process is obscure, this choice 
is arbitrary, and we have set H to be independent of all 
parameters. The relation between the terms is iJo(Ferland) 
= i/(Bayet)/(n/l cm~^), or, equivalently, -ffo(Ferland) = 
if (Bayet)(T/4 x 10^ K cm^'^) for the assumed pressure. 
Cosmic ray induced ionisation of a region with a hydrogen 
number density of 10"^ cm^"^ will be subject to a comparable 
heating rate (i.e. H = lO^^^erg cm^'^ s^^) due to cosmic 
rays if C is within an order of magnitude of 10^^^ s^^. 



Table 1. Summary of m.odcl runs. When the symbol for a molecular species is underlined, its fractional abundance is above 1x10^^" otherwise it is simply above 
1x10^^^ (limit of detectability assumed). The temperature mentioned in the column 6 is self-consistently determined at the particular value of Av shown. 

Model 



coTh 

1.7x10^'' 
4.1x10-" 
6.5x10"" 
3.5x10"" 
1.2x10"^ 

9.9x10"" 
1.0x10"^° 
3.9x10"" 
3.8x10"" 
5.35x10"" 
5.8x10"" 
3.7x10"" 
3.7x10"" 
1.4x10"^ 

1.4X10"-''' 

1.3x10"'' 

l.SxlO""' 

1.3x10"'' 



Av = 3 

1 
2 
3 

4 
5 

6 

7 
8 
9 
10 
11 
12 
13 
14 

15 

16 

17 

18 



H 

(erg cm"'' s"^) 

1x10"^" 
1x10"^^ 
1x10"^" 
1x10"^^ 
1x10"^^ 

1x10"^^ 
1x10"^^ 
1x10"^° 
1x10"^^ 
lxlO"^-°'-lxlO"" 
1x10"^" 
1x10"^^ 


1x10"^^ 



1x10"^" 
1x10"^^ 
1x10"^" 



P 

(cm"^ K) 

4x10" 
4x10" 
4x10" 
4x10" 
4x10" 



I c 

(Habing) (s-i) 



100 
100 
100 
100 
100 



5x10"" 
5x10"" 
5x10"^'' 
5xl0"^-°' 
5x10"^^ 



T 
(K) 

2866.9 
91.4 
68.2 
51.3 
25.3 



xlO" 


10 


5x10" 


13 


2479.6 


xlO" 


10 


5x10- 


13 


2456.1 


xlO" 


10 


5x10" 


14 


87.5 


xlO" 


10 


5x10" 


14 


86.2 


xlO" 


10 


5x10" 


15 


10000.0 


xlO" 


10 


5x10" 


15 


61.8 


xlO" 


10 


5x10" 


15 


51.0 


XlO" 


10 


5x10" 


15 


50.8 


xlO" 


10 


5x10" 


17 


13.8 



4x10" 
4x10' 
4x10' 
4x10" 



10 
10 
10 
10 



5x10"'' 
5xl0"^= 
5 XlO"'-"' 
5x10"^' 



12.4 
43.8 
41.1 

74.8 



Interesting 
molecules 

O, OH= H3+, C;!:, CN, C, H2, H, H2O, CO 

H3+ , C+ , CN , C, H2, H, H2O , HCN, C2, H3O+ , C2H, HCO+ . CO , O, OH, CO2, HNC 

H3+, Ci, CN, C, H2, H, H2O , HCN, NH3, C2, H3O+ , C2H , OCN, SO, CS, HCO+ , CO , O, 

hI^ , C+ , CN, C, H2, H, H2O , HCN, NH3, C2, H3O+ , C2H, OCN, SO, CS, HCO+ , CO, O, 

H3+ , C+ , CN, HCO, C, H2, H, H2O. HCN , NH3, C2., H3O+, C2H, OCN, HCS, H2CS, NS, 

O, OH, H2CO, HS, CO2, HNC . HNO 

H3+ , C+ , C, H2, H, H2O, CO, O, OH 

H3+ , C+ , C, H2, H, H2O, CO, O, OH 

C+ , CN, C, H2, H, H2O, HCN, C2, H3O+ , C2H, HCO+ , CO, O, OH, CO2, HNC 

H3+ , C+ , CN, C, H2, H, H2O, HCN, C2, H3O+ , C2H, HCO+ . CO, O, OH, CO2, HNC 

Ql, C, H2, H, O, OH 

H3+ , C+ , CN, C, H2, H, H2O, HCN, NH3, C2, H3O 

H3+ , C+ , CN, C, H2, H, H2O , HCN, NH3 

H3+ , C+ , CN, C, H2, H, H2O , HCN, NH3 

H3+ , C+ , CN, HCO, C, H2, H, H2O, HCN , NH3 , C: 
HCO+, CO, O, OH, H2CO, HS, CO2 , HNC , HNO 
H3+ , C+ , CN, HCO, C, H2, H, H2O, HCN , NH3 , C: 
HCO+, CO, O, OH, H2CO, HS, CO2 , HNC , HNO 
H3+ , C+ , CN, C, H2, H, H2O, HCN, NH3, C2, H2S, H3O 
HCO+ , CO , O, OH, H2CO , HS, CQ 

H3+ , C+ , CN, C, H2, H 

HCO+ , CO, O, OH, H2CO , HS, CO2 , HNC , HNO . C2N, H2CN 

H3+ , C+ , CN, C, H2, H, H2O, HCN, NH3, C2, H3O+, C2H, OCS, OCN, HCS, H2CS, NS, 

HCO+, CO. O. OH, H2CO, HS, CO2, HNC, HNO, H2CN 



, OH, H2CO, CO2, HNC . HNO m 
OH . H2CO, CO2, HNC, HNO ^ 
SO, CS, HCO+, CO, ^• 



o 



HaO^ 
HaO^ 



C2H, OCN, SO, CS, HCO+ , CO . O, 
C2H, OCN, SO, CS, HCO+ , CO, O, 
C2H, OCN, SO, CS, HCO+ , CO, O, 
H2S, H3O+, C2H, OCS, OCN, HCS, 



OH . H2CO, CO2, HNC . 
OH . H2CO, CO2, HNC, 
OH . H2CO, CO2, HNC, 
H2CS, NS, SO, CS, 



HNO 
HNO 
HNO 



H2S, H3O+, C2H, OCS, OCN, HCS, H2CS , NS, SO, CS, 
OCS, OCN. HCS, H2CS. NS. SO, CS. 



C2H 

HNC , HNO , C2N, H2CN 
H2O, HCN, NH3, C2, H2S, H3O+, C2H, OCS, OCN. HCS, H2CS 



, NS, SO, CS, 
SO, CS, 



I 
o 






Av = 8 

8 
9 

10 
11 
12 
13 
14 

15 

16 

17 

18 



1x10"^" 
1x10"^^ 
lxlO"^-°'-lxlO" 
1x10"^° 
1x10"^^ 


1x10"^^ 



1x10"^" 
1x10"^^ 
1x10"^° 



4x10" 
4x10" 
4x10" 
4x10" 
4x10" 
4x10" 
4x10" 

4x10" 

4x10' 

4x10' 

4x10" 



10 
10 
10 
10 
10 
10 
10 

10 

10 

10 

10 



5x10"" 
5x10"" 
5x10"^'' 
5x10"^'' 
5x10"^"' 
5x10"^"' 
5x10"^' 

5x10"" 

5xl0"^= 

5x10"^"' 

5x10"" 



101.7 
99.1 
10000.0 
78.7 
60.2 
59.9 
15.0 

13.4 

47.5 

45.2 

63.3 



4.9x10"" 
4.7x10"" 
5.36x10"" 
8.3x10"" 
5.6x10"" 
5.5x10"" 
1.4x10"-^ 

1.4x10"'' 

1.3x10"'' 

1.3x10"'' 

1.4x10"'' 



H3+ 


C+, CN, 


H3+ 


C+, CN, 


c+. 


C, H2, H, 


H3+ 


C^CN, 


H3 + 


C+, CN, 


H3 + 


C+, CN, 


H3 + 


C+, CN, 


HOC 


)+, CO, 


H3+ 


C+, CN, 


HCC 


)+, CO, 


H3+ 


C+, CN, 


HCC 


)+, CO, 


H3+ 


C+, CN, 


HCC 


)+, CO, 


H3+ 


C+, CN, 


HCC 


)+, CO, 



C_, H2, H, H; 
C- , ri2 , £1, M; 

O, OH ~ 

C- , ri2 , £1, M; 

C , H2, H, H: 

C- , ri2 , £1, M; 

HCCX C, H^, 
, OH, H2CO, 
HCO, C, H2, 
, OH, H2CO, 
C, H2, H, H2 
, OH, H2CO , 
C, H2, H, H; 
, OH, H2CO, 
C, H2, H, H2 
, OH, H2CO 



O, HCN, 
O, HCN, 

O, HCN, 
O, HCN, 
O, HCN, 
H, H2O. 
HS, CO2 
H, H2O. 
HS, CO2 
O, HCN , 
HS, CO2 
O, HCN , 
HS, CO2 
O, HCN , 
HS, CO2 



NH; 
NH< 



, C2 
, C2 



H30^ 
H30^ 



, C2H, SO, HCO+ 
, C2H, SO, HCO+ 



CO . O. OH, CO2, HNC 
CO . O, OH, CO2, HNC 



a. 



, CO, O, OH, H2CO, CO2, HNC . 
, CO, O, OH, H2CO, CO2, HNC . 
C2H, OCN, SO, CS, HCO+ , CO, O, OH, H2CO, CO2, HNC . 
H3O+, C2H, OCS . OCN, H2CS, CH3OH, NS, SO, CS, 

H2S, H3O+, C2H, OCS . OCN, H2CS, CH3OH, NS, SO, CS, 

C2H, OCS, OCN . HCS, H2CS, NS, SO, CS, 



NH3, C2, H3O+ , C2H, OCN, SO, CS, HCO 

NH3, C2, H3O+ , C2H, OCN, SO, CS, HCO 

NH3, C2, H3O+ 

HCN , NH3 , C2, H2S 

, HNC , HNO 

HCN , NH3 , C: 

, HNC , HNO 

NH3, C2, H2S, H3O 

, HNC , HNO , C2N 

NH3, C2, H2S, H3O 



HNO 
HNO 
HNO 



C2H, OCS, OCN. HCS, H2CS, NS, SO, CS, 



, HNC , HNO , C2N 
NH3, C2, H2S, H3O+ 
, HNC , HNO 



C2H, OCS, OCN .H2CS, NS, SO, CS, 



13 
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The higher values of our assumed heating rates per unit 
volume and cosmic ray induced ioni sation r ates roughly 
approach the values found by Ferla nd et alj ([2008, 2009) 
to be compatible with the optical emission line spectra in 
NGC 1275. However, we limited the upper range of the as- 
sumed rates to avoid heating the gas beyond the maximum 
temperatures at which molecules are abundant. This limit- 
ing is consistent with the values of H and C, decreasing with 
depth into a filament. 

3. Results 

The code determines self-consistently the steady-state ther- 
mal and chemical properties of a constant-pressure fila- 
ment subjected to the three heating processes described in 
Section[21 of these, H and C affect the bulk of the filament. 
In contrast, the effects of / are confined to narrow surface 
regions, as a comparison of results given in Figures [5] and [3] 
for our models 14 and 5 (which differ only in /) confirms. 
The thermal and chemical properties are computed as func- 
tions of Av, measured from the outer edge of the filament. 
We present in Table [1] a summary of the models explored 
and the corresponding thermal and chemical properties at 
two values of Av, 3 and 8 mag. These values are intended to 
represent conditions near the edge and in the dark interior 
of a filament. We do not list results at Av = 8 for Models 
1-5 as they should not differ from those at Av = 3. We in- 
clude in the final column of Table 1 a list of molecules that 
are predicted to be fairly abundant in each model (under- 
lined species) and molecules that may be near the current 
limit of detection (non- underlined species) . Evidently, some 
models are expected to be chemically very rich, while oth- 
ers should be chemically poor. Table [2] contains values at 
Av = 8 of the fractional abundances of the species. In this 
section we focus on the results shown in various figures, but 
the trends that we infer from those figures are also seen in 
the results in Tabled] 

The chemistry varies little with Av as Av increases be- 
yond 8. However, the figures that we show give results to 
higher values of Av because the visual extinctions of the 
clumps of which a filament is composed may range to higher 
values. The maximum Av from its edge to its centre is prob- 
ably that of the maximum mass isothermal sphere that is 
stable against self-gravity. For clumps with dust and ele- 
mental abundances like those of of clouds in the local in- 
terstellar medium, this Av is given roughly by 



Av^a. = 15mag[n(H2)r/(4 x lO^Kcm"^)]^/^ 



(1) 



This result is based on the expression for the maximum 
external pressure that can be applied to a Bonner-Ebert 
sphere without causing it to collapse, the assumption that 
the pressure in the sphere is roughly constant, and that 
a column density of 2 x 10^^ cm~^ hydrogen nuclei corre- 
sponds to one magnitude of visual extinction. 

The temperature quoted in Table [T] is the kinetic tem- 
perature determined self-consistently with the chemistry, 
given for depths into the filament corresponding to these 
two values of Av. The fractional abundance of CO deter- 
mined by the code at those same depths is also given, as a 
crude measure of the effectiveness of the chemistry in the 
filament. It ranges from a value of 1.3 x 10~^, correspond- 
ing to a rich chemistry in which about 10% of available 
carbon is in CO, to very small values. We summarise in 



Table [T] the chemistry determined for each model by sim- 
ply listing molecules familiar from studies of galactic inter- 
stellar media; more species than those listed are present. 
Molecules listed have a fractional abundances of at least 
1 X 10~^^ (adopted as an arbitrary limit of detectability), 
and those underlined have fractional abundances of at least 
1 X 10"^". Detections of C2 and CO2 will be unlikely, but 
these molecules are listed for completeness. More informa- 
tion about the dependence of gas density, temperature and 
chemistry on Av for particular models is found in Figures 
[U [5] and m described below. 

Variations in the additional heating rate: 

Model 10 explores the effect of varying the additional 
rate, H , over the range 1 x 10~^^ to 1 x 10~^^ erg cm~^ s~^. 
The high values of H dominate other heating mechanisms 
and overwhelm the cooling mechanisms. The resulting tem- 
perature is so high that very few molecules can survive. 
Since H is assumed to be independent of depth, this result 
is true at all depths. 

For values of iJ of 1 x 10^^° erg cm"'^ s~^ and smaller, 
a rich chemistry can exist but is sensitive to the adopted 
cosmic ray ionisation rate. For example. Models 11 {H — 
1 x 10^^° erg cm~^ s~^) and 12 (iJ = 1 x 10^^^ erg cm^^ 
s~^) are seen in Tabled] to have almost identical extensive 
chemistries, while Model 13 [H — 0) is indistinguishable 
from Model 12. Thus, H = Ix 10^^° erg cm~^ s~^ appears 
to be a critical value; above this value, the heating is so 
great that the chemistry is largely suppressed, while below 
this value a rich chemistry occurs and is insensitive to H. 
However, as mentioned above, the value of ( affects the 
richness of the chemistry for a particular value of H. For 
example. Models 11 and 8 (both with H = 1 x 10~^° erg 
cm-3 s-i) differ in C (5 x 10"^^ s^^ and 5 x 10"^'' s"\ 
respectively); Model 11 is distinctly richer in chemistry. 

For heating due to dissipation to exceed that due to 
cosmic ray induced ionisation H must exceed a value given 
roughly by 



H„,^n = lxl0"^Vgcm"3s"ix(n(H2)T)/(4xl0^Kcm"3)x(2) 



(T/WOK)-^ X C/(5 X 10""s"i). 

This result is based on the assumption that one ionisa- 
tion produces about 40 eV of thermal energy. 

Figure [1] shows the dependences of the total number 
density and temperature in the gas, and the chemical abun- 
dances, as functions of depth for Model 16. This model may 
have plausible values of the physical parameters. It shows 
a modest temperature and a rich chemistry. Evidently, the 
familar galactic tracers CS, SO, HCO+, C2H, CN, HNC 
and HCN should be useful for determining the nature of 
the intracluster material in cases where the parameters are 
similar to those of Model 16. 

Figure [2] compares the chemistries predicted by two 
Models, 14 and 18, which both have low values of C but 
differ strongly in 7J (1 x 10~^^ erg cm~'^ s~^ and 1 x 10"^" 
erg cm~^ s~^, respectively). Some spikes are seen for Model 
18 because the heating is only just matched by the cool- 
ing, and the balance is thus computationally delicate. Both 
models show a rich chemistry, but there are some physi- 
cal and chemical differences. For example, the important 
tracer ion HCO+ is more abundant in Model 18; and the 
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Fig. 1. Fractional abundances with respect to total hydrogen number density of various species (X) as functions of optical 
depth (Ay), derived from model 16. For each model presented in Figs [TJ [5] and [31 the plots on the top correspond to 
the density and the temperature profiles. The chemical formula of each molecular species studied is indicated on the 
plots. When not mentioned, it means that the fractional abundance of X is lovifcr than the arbitrary limit of detectability 
assumed to be equal to 1 x 10""'^^. We have used solid lines for representing the model having an additional heating source 
of 1 X 10^^° erg cm^'^ s~^ and the dashed lines show the models having an additional heating source of 1 x 10^^^ erg 
cm ■^ s . 



molecules CN, HCN and HNC have a different priority or- 
der of abundance in the two models, caused by the higher 
temperature in Model 18. 

Variations in the cosmic ray ionisation rate: 

The value of C is an important parameter both for 
the chemistry and for the thermal balance. The cases 
with highest C examined here are Models 6 and 7, with 



( — 5 X 10~^^s~^. The value of H is low in both Models 
(1 X 10~^^ erg cm~^ s~^ and 1 x 10"^'^ erg cm~^ s~^, re- 
spectively), so the major heat input is from the cosmic rays. 
The temperature is high and the chemistry is largely sup- 
pressed; even CO is held down to a fractional abundance 
of about 1 X 10^^°. The high C but low CO abundance 
leads to a significant Hg abundance. This suppression of 
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Fig. 2. Fractional abundances with respect to total hydrof 
depth (Av), derived from Model 18 (left) and Model 14 
having an additional heating source of 1 x 10^^° erg cm^ 
the models having an additional heating source of 1 x 10 



jen number density of various species (X) as functions of optical 
(right). We have used solid lines for representing the models 
^ s~^ (i.e Model 18) and, as previously, the dashed lines shovif 
-22 gj.g (.j^-3 g-i j^jg Model 14). See caption of FigH] 



th e chemistry for high C is in agree ment with the findi ngs 
of lLepp fc Dalgarncl (jlQQeD fsee also lLepp fc Tinelll998[ ). 

Figure [3] compares results from Models 5 and 12 which 
both have low values of i/ (1 x 10"^^ erg cm"-^ s^^), but 
strongly differing values of C (5 x 10~^^ s~^ and 5 x 10~^^ 
s~^, respectively). Both are chemically rich (Model 5 is 
sHghtly richer) but there are significant differences. For ex- 
ample, Model 5 has very much higher SO and CS abun- 
dances than Model 12, while the reversed is true for C2H 
and HCO"'". The abundance order of priority of CN, HNC 
and HCN is reverse between Models 5 and 12, with the 



abundance of CN being one order of magnitude larger in the 
high ^ case. Model 12. This case also has abundant atomic 
species H, C and C"*" relative to Model 5. It is well known 
that the H:H2 ratio, maintained in this c ase by the high 
C, has a significant effect on the chemistry (iRawlines et al.l 
2002). 

Other variations: 

The adopted pressure, P, has a limited effect on the 
thermal and chemical behaviour. Comparing Models 11 and 
16, or 12 and 17 (which differ only in P) shows that the 
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Fig. 3. Fractional abundances with respect to total hydrogen number density of various species (X) as functions of optical 
depth (Av), derived from Model 12 (left) and Model 5 (right). As previously, the dashed lines show the models having an 
additional heating source of 1 x 10^^^ erg cm^"^ s^^ (i.e both Models 12 and 5). Here the difference between the models 
is the cosmic ray ionisation rate used (^ = 5 x 10~^^ s~^ for Model 12 whereas Model 5 has ( = 5 x 10"^'' s""'^) and the 
FUV radiation field (/ = 10 Habing for Model 12 whereas Model 5 shows / = 100 Habing). See caption of FigUand text 
in Sect. [31 



temperature is slightly lower at higher pressure (as expected 
since the cooling is more effective), and that this slight dif- 
ference generates a difference in the chemistry, with the 
cooler models being very slightly chemically richer. 

The ambient radiation field can be an important source 
of heat at the cloud edge, comparable to if = 1 x 10~^° 
erg cm~^ s~^ in unshielded situations. However, this heat- 
ing effect decays rapidly with depth into the cloud. For 
the higher values of / explored here, / = 100, the radia- 



tion field may still be making a contribution to the heat- 
ing, even at Av = 3. We show in Table |3] a comparison 
of the heating rates from the main heating sources in sev- 
eral models at both low and high Av. Evidently, the cosmic 
rays and the dissipation (i.e. H) are the most important 
sources. Comparing Models 4 and 12 shows that the larger 
value of / in Model 4 can be making no difference when 



C = 5x 10- 



-15 



because cosmic rays dominate the heat- 



ing in these cases. However, comparing Models 5 and 14, 
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Table 3. Main contributors to the heating rates (erg s ^ 
cm~^) in several models (See Table [Ij. 



Av = 3 mag- 
heating sources 
FUV 
H 
Cosmic rays 



Model 1 
8.8 X 10"^* 
1.0 X 10^^° 
1.0 X 10"^^ 



Model 8 
8.8 X 10"^^ 
1.0 X 10"^" 
3.4 X 10"^" 



-23 



Model 18 
1.7 X 10 
1.0 X 10" 
4.0 X 10" 



Av = 8 mag 

heating sources 

FUV 

H 

Cosmic rays 



Model 14 Model 8 

8.8 X 10"^® 7.6 X 10"^^ 
1.0 X 10"^^ 1.0 X 10"^° 



2.0 X 10" 



2.9 X 10" 



Model 18 
1.3 X lO"^'^ 
1.0 X 10"^° 
4.7 X 10"^=* 



both of which have a low value of C, 5 x 10~^^ s~^, we see 
that the higher value of / in Model 5 increases the temper- 
ature relative to Model 12 - though not enough to make a 
difference to the chemistry. 

At low Av, Model 1 shows a surprisingly high temper- 
ature and poor chemistry compared to a similar model, 8, 
with lower value of /. Evidently, the radiation field can 
have a significant effect in the boundary layers of filaments. 
This comparison emphasises the "fiip-flop" effect; a rich 
chemistry provides many coolants that generate a low tem- 
perature, enhancing the chemistry; and vice-versa. Model 
1 is expected to show a much lower temperature at greater 
depths. 



4. Conclusions 

The self-consistent treatment of thermal and chemical be- 
haviour tends to lead to one of two alternative regimes: low 
temperature and rich chemistry, or high temperature and 
poor chemistry. The rich chemistry provides efficient cool- 
ing, enhancing the chemistry, while elevated temperatures 
suppress the chemistry and destroy coolants. 

The parameter regime that is covered in this study 
shows that there are two main contributors to the heat- 
ing, the additional source H and cosmic rays. The ambient 
radiation field is unlikely to be dominant and is in any case 
confined to a skin-effect of a few visual magnitudes into 
the filament. The density of the filament has been assumed 
in our calculations to be uniform. Of course, is the gas is 
clumpy, then the radiation fields may play a larger role at 
greater depths than in the homogeneous case. Varying the 
pressure from the adopted value has only a minor effect on 
the chemistry. 

Very high values oi H , > 1 x 10^^° erg cm^^ s^^, tend 
to suppress the chemistry and consequent coolants, so tem- 
peratures approach those of HII regions. Very high cosmic 
ray ionisation rates, > 1 x 10"^'^ s~^, have a similar ef- 
fect. For iJ < 1 X 10-20 erg cm'^ s'^ and C < 1 x IQ-^^ 
s""'^, the heating and consequent chemical effects of the two 
processes combine. 

These rates considerably exceed the typical rates in dark 
clouds near to the Sun (and may be even higher in the fila- 
mentary optical regions) . If these are representative of the 
conditions in filamentary molecular regions, then there are 
many molecular species that should be abundant and use- 
ful as tracers of the physical conditions. CS, CN and HNC 
decline less in abundance with increasing heating and ioni- 
sation rates than many species. C2H and HCO"'' increase in 
abundance with cosmic ray ionisation rate for a wide range 



of rates but become less abundant for very high rates. The 
ratio of those two abundances decreases with increasing dis- 
sipative heating rate. N2H''", a molecule used successfully 
to trace dense molecular clouds in the Galaxy, is notable 
by its absence from our predictions. 

In our recommendations for molecular line searches in 
cluster filaments, we have been guided by the molecular 
abundances predicted by our computations and by their 
sensitivity to changes in heating and ionisation rates. The 
construction of detailed models will be possible once more 
observational data exist. However, we can give some indi- 
cation of the expected range of HCO"'", C2H and CN line 
intensities, assuming no bea m dilution. To do so, we have 
used RADEXQ developed bv Ivan der Tak et all ()2OO70 . for 
plane-parallel geometry, for Models 12, 14, 16 and 18 seen in 
Figs. El and [31 We have obtained for the HCO+(1-0) tran- 
sition, values ranging from 0.62 K knis~^ (Model 14) to 
5.6 K kms-i (Model 16). For the HCO+(5-4) line, values 
range from 0.2 K kms"^ (Models 12 and 18) to 4.2 K kms"^ 
(Model 16). For the CN(l-O) line, we have obtained values 
from 0.2 K kms'^ (Model 14) to 7.6 K kms"^ (Model 16). 
For the CN(4-3) lines, we have obtained values ranging from 
1.4x10^2 K kms-i (Model 18) to 1.6 K kms"! (Model 16). 
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Table 2. Fractional abundances with respect to the total number of H-nuclei for models 8-18 at Av = 8 mag. The power 
is indicated in parenthesis such as a(b)=axlO^. When a dash is used, it means that the fractional abundance is below 
our assumed limit of detectability of 1x10" 
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Model: 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


Species 
























Hj 


2.84(-8) 


2.76(-8) 


- 


5.36(-8) 


3.19(-8) 


3.16(-8) 


9.09(-10) 


8.35(-10) 


1.34(-8) 


1.27(-8) 


2.74(-9) 


C+ 


6. 11 (-6) 


6.18(-6) 


1.40(-5) 


1.70(-6) 


1.96(-6) 


1.97(-6) 


6.01(-10) 


5.32(-10) 


7.58(-8) 


7.43(-8) 


3.27(-9) 


CN 


9.56(-ll) 


9.33(-ll) 


- 


2.28(-9) 


1.99(-9) 


1.99(-9) 


1.81(-10) 


1.53(-10) 


7.24(-9) 


7.23(-9) 


8.01(-10) 


HCO 


- 


- 


- 


- 


- 


- 


2.61(-12) 


2.56(-12) 


- 


- 


- 


C 


3.23(-6) 


3.36(-6) 


2.06(-7) 


4.19(-6) 


6.66(-6) 


6.73(-6) 


4.48(-9) 


4.05(-9) 


1.07(-6) 


1.08(-6) 


3.15(-8) 


H2O 


2.30(-7) 


2.18(-7) 


- 


9.90(-8) 


5.29(-8) 


5.22(-8) 


4.97(-7) 


5.14(-7) 


1.81(-7) 


1.78(-7) 


3.82(-7) 


HON 


2.87(-12) 


2.68(-12) 


- 


9.25(-ll) 


4.58(-ll) 


4.51(-11) 


1.60(-9) 


2.02(-9) 


1.21(-9) 


1.20(-9) 


4.61(-10) 


NH3 


1.22(-12) 


1.13(-12) 


- 


1.15(-10) 


4.50(-ll) 


4.42(-ll) 


5.46(-8) 


6.85(-8) 


3.12(-9) 


3.02(-9) 


1.66(-8) 


H2S 


- 


- 


- 


- 


- 


- 


3.70(-ll) 


4.07(-ll) 


1.54(-12) 


1.60(-12) 


l.Ol(-ll) 


H3O+ 


3.14(-9) 


2.94(-9) 


- 


8.10(-10) 


4.23(-10) 


4.16(-10) 


2.41(-11) 


2.33(-ll) 


1.89(-10) 


1.77(-10) 


5.11(-11) 


C2H 


1.98(-12) 


2.00(-12) 


- 


2.84(-10) 


2.95(-10) 


2.95(-10) 


5.21(-12) 


5.11(-12) 


2.74(-10) 


2.76(-10) 


3.43(-12) 


ocs 


- 


- 


- 


- 


- 


- 


1.21(-10) 


1.50(-10) 


2.05(-ll) 


2.19(-11) 


1.125(-11) 


OCN 


- 


- 


- 


8.66(-12) 


3.80(-12) 


3.74(-12) 


2.54(-8) 


3.26(-8) 


4.82(-10) 


5.07(-10) 


3.52(-10) 


HCS 


- 


- 


- 


- 


- 


- 


- 


- 


3.55(-ll) 


3.74(-ll) 


- 


H2CS 


- 


- 


- 


- 


- 


- 


6.82(-ll) 


6.72(-ll) 


2.01(-10) 


2.13(-10) 


4.92(-ll) 


CH3OH 


- 


- 


- 


- 


- 


- 


1.99(-12) 


2.44(-12) 


- 


- 


- 


NS 


- 


- 


- 


- 


- 


- 


3.54(-12) 


2.96(-12) 


5.64(-12) 


5.42(-12) 


1.13(-11) 


SO 


1.21(-12) 


1.08(-12) 


- 


2.03(-ll) 


4.84(-12) 


4.70(-12) 


3.60(-8) 


3.56(-8) 


1.14(-9) 


1.085(-9) 


2.03(-8) 


cs 


- 


- 


- 


3.65(-10) 


3.59(-10) 


3.58(-10) 


4.76(-8) 


4.67(-8) 


1.86(-8) 


1.93(-8) 


4.64(-8) 


HCO+ 


4.05(-10) 


3.72(-10) 


- 


8.80(-10) 


2.94(-10) 


2.87(-10) 


5.32(-ll) 


4.75(-ll) 


5.94(-10) 


5.39(-10) 


2.62(-10) 


CO 


4.86(-6) 


4.66(-6) 


- 


8.31(-6) 


5.57(-6) 


5.50(-6) 


1.39(-5) 


1.38(-5) 


1.30(-5) 


1.30(-5) 


1.41(-5) 





2.53(-5) 


2.56(-5) 


3.18(-5) 


2.29(-5) 


2.60(-5) 


2.60(-5) 


9.13(-6) 


9.25(-6) 


1.73(-5) 


1.73(-5) 


1.18(-5) 


OH 


1.48(-6) 


1.40(-6) 


1.67(-12) 


4.08(-7) 


2.23(-7) 


2.20(-7) 


2.58(-8) 


2.59(-8) 


1.17(-7) 


1.10(-7) 


3.77(-8) 


H2CO 


- 


- 


- 


2.93(-ll) 


2.23(-ll) 


2.21(-11) 


5.13(-11) 


5.06(-ll) 


5.18(-10) 


5.25(-10) 


3.92(-ll) 


HS 


- 


- 


- 


- 


- 


- 


1.23(-11) 


1.21(-11) 


8.89(-12) 


8.78(-12) 


1.82(-11) 


HNC 


1.22(-11) 


1.15(-11) 


- 


2.27(-10) 


1.31(-10) 


1.29(-10) 


9.33(-10) 


1.15(-9) 


1.41(-9) 


1.40(-9) 


4.22(-10) 


HNO 


- 


- 


- 


6.61(-11) 


4.93(-ll) 


4.88(-ll) 


4.22(-10) 


2.84(-10) 


1.99(-9) 


1.95(-9) 


7.43(-10) 


C2N 


- 


- 


- 


- 


- 


- 


- 


- 


1.77(-12) 


1.75(-12) 


- 
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